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We show that relativistic fluids behave as non-Newtonian fluids. First, we dis- 
cuss the problem of acausal propagation in the diffusion equation and introduce 
the modified Maxwell-Cattaneo-Vernotte (MCV) equation. By using the modified 
MCV equation, we obtain the causal dissipative relativistic (CDR) fluid dynamics, 
where unphysical propagation with infinite velocity does not exist. We further show 
that the problems of the violation of causality and instability are intimately related, 
and the relativistic Navier-Stokes equation is inadequate as the theory of relativistic 
fluids. Finally, the new microscopic formula to calculate the transport coefficients 
of the CDR fluid dynamics is discussed. The result of the microscopic formula is 
consistent with that of the Boltzmann equation, i.e., Grad's moment method. 



I. INTRODUCTION 



Typical examples of fluid are water and air, whose dynamics is described by the Navier- 
Stokes (NS) equation. These are called Newtonian fluids. There are, however, various fluids 
called non-Newtonian fluids, which cannot be described by the NS equation. The difference 
between these two types of fluid comes from the behavior of the shear stress tensor. In Fig. 
[TJ the various shear stress tensors are shown as a function of the gradient of the fluid velocity. 
When the shear stress tensor increases proportionally with the velocity gradient, the fluid 
is Newtonian, which is denoted by the solid line. Non-Newtonian fluids exhibit a more 
complex behavior as is shown by the dashed lines. The Bingham flow 1) represents a similar 
linear relation but the shear stress tensor does not disappear even in the vanishing velocity- 
gradient limit. The dilatant fluid 2) and pseudoplastic 3) show non-linear dependences. The 
shear stress tensor of the thixotropic fluids 4) depends on time. 

If the dynamics of relativistic many-body systems can be described by using coarse- 
grained equations such as fluid dynamics, is the behavior of relativistic fluids Newtonian or 
non-Newtonian ? In order to answer this question, we will start our discussion from diffusion 
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FIG. 1: Shear stresses of Newtonian and non-Newtonian fluids as a function of the velocity gradient. 
The solid line denote a Newtonian fluid. The dashed lines represents 1) Bingham flow, 2) dilatant 
fluids, 2) pseudoplastic and 4) thixotropic fluids, respectively. 

processes, because the problem which we will encounter in relativistic fluid dynamics has 
already appeared in the diffusion equation. 



II. DIFFUSION EQUATION AND MCV EQUATION 

We consider a random walk process, where a particle moves to left or right on a one- 
dimensional lattice with equal probability. The probability distribution function P(x, t) 
satisfies P(x,t + dt) = (P(x — dx, t) + P(x + dx, t))/2, where dx and dt are the size of lattice 
and time steps, respectively. In the continuum limit, we obtain the diffusion equation, 

dtP(x,t) = Dd 2 x P(x,t), (1) 

with the definition of the diffusion coefficient, D = \im.dt,dx-+o(dx) 2 /dt. Then the particle can 
move by dx at each time step dt and the velocity of the particle is given by v — dx/dt. Note 
that the continuum limit should be taken by fixing D. This leads to the infinite velocity, 

lim v = lim y^D/dt = oo. (2) 

dt,dx-^0 dt,dx^0 

Let us consider a possible modification of the diffusion equation to avoid this violation 
of causality. Remember that the diffusion equation consists of two structures. One is the 
equation of continuity, 

d t n + VJ = 0, (3) 
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where n is a conserved density and J is a current. The other is the definition of J. To obtain 
the diffusion equation, we assume that J is proportional to the corresponding thermodynamic 
force F, 

J = —DF = -DVn. (4) 

Here F = Vn for the diffusion process. This is called Fick's law. 

The equation of continuity ^ should be satisfied for any conserved density. Thus, if it 
is possible to derive a modified diffusion equation consistent with causality, only Eq. Q 
can be changed. As a matter of fact, from a microscopic theory such as the linear response 
theory, a more general expression of J is given by the time convolution integral, 

J(t) = - J dsG(t- s)F(s), (5) 

where G(t) is the memory function which is given by the time correlation function of mi- 
croscopic degrees of freedom. Thus the time scale of the memory function is characterized 
by the microscopic time scale. If the time scale of macroscopic variables such as J and F is 
clearly separated from the microscopic one, we can approximately replace the time depen- 
dence of G(t) with the Dirac delta function, G(t) = DS(t) and then we can reproduce Fick's 
law Q. 

When, however, the time scales are not clearly separated, the time dependence of G(t) 
should be taken into account. As a simplest choice, we use the exponential form, 

G(t) = -e-*'**, (6) 

where tr is the relaxation time which characterizes the microscopic time scale. Substituting 
into Eq. ^ and operating the time derivative, we obtain 

T R d t 3(t) + 3(t) = -DF(t). (7) 

This is the so-called Maxwell-Cattaneo-Vernotte (MCV) equation. When there is a clear 
separation of microscopic and macroscopic time scales, tr vanishes and the MCV equation 
is reduced to Fick's law Q. 

By eliminating J from the equation of continuity ^ with the MCV equation ([7]), we have 
the telegraph equation for n, 

r R d 2 t n + d t n = DV 2 n. (8) 
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FIG. 2: The time evolutions of the diffusion equation (dashed line) and the MCV equation (solid 
line) for t/rn = 0.5, 2 and 6. The initial condition is shown by the dotted line at £/tr = 0.5. 

The analytic solution is given in Ref. [Ij. For example, the solution of 1+1 dimensional 
system is given by 

e -t/(2r R ) 

n(x,t) = {n (x + vt) + n (x — vt)} 



+e 



+e 



"X+Vt 

- !i2 ^ ] I dx' 

x—vt 
x+vt 



_L no(x ') + h x J {x')\ J {y/vH*-{x- X ')*/{2vT R )) 



*/( 2t «) / dx' cVo ( \/v 2 t 2 - (x - x') 2 /(2vt r ) 



(9) 



l x - vt 2vt r 

where Iq is the modified Bessel function, and n and Jo are initial conditions of the conserved 
density and corresponding current, respectively. The maximum propagation speed of this 
equation is characterized by v, 

v = y/Dfa, (10) 
which diverges in the diffusion limit (tr = 0). In Fig. [2j the time evolutions of the diffusion 
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equation and the telegraph equation are shown by the dashed and solid lines. We use 
D — tr — 1, leading to v — 1. The initial denstiy distribution n (x) is plotted in the figure 
of t/rR = 0.5 by the dotted line, and we further set Jo(x) = 0. One can see that, because of 
the memory effect, there exists a non-trivial structure at the boundary of the expanding n, 
reflecting the initial distribution in the telegraph equation, although the diffusion equation 
always shows the Gaussian forms. Of course, because of the finite propagation speed, the 
expansion of the telegraph equation is slower than that of the diffusion equation. See also 
Ref. 0. 

In table [T| the comparison of the diffusion equation and the MCV equation is summarized. 
If the diffusion equation is a coarse-grained dynamics of the underlying microscopic physics, 
it should be derived from a microscopic theory using systematic coarse-graining. As a 
matter of fact, as is discussed in textbooks, it is believed that the diffusion equation can be 
derived with the projection operator method. However, we should note that one non-trivial 
approximation is used in this derivation. As a matter of fact, it was recently found that, 
when such non-trivial approximation is not applied, the MCV equation is obtained instead 
of the diffusion equation [3]. 

Correspondingly, the microscopic Hamiltonian which describes a diffusion process should 
have a symmetry associated with the conserved density, and we can derive the corresponding 
sum rule. The sum rule determines the initial time evolution of the conserved density. The 
telegraph equation is consistent with this sum rule, although the diffusion equation is not 

As is well known in linear irreversible thermodynamics, the positivity of the entropy 
production is algebraically satisfied when there is a simple linear relation between J and 
F. However, this linear relation is not satisfied in the MCV equation. Moreover, when 
there is no clear separation of time scales, we cannot assume quasi-adiabatic changes of 
thermodynamical variables, and then heat will play a more fundamental role instead of 
entropy. Thus the second law for the MCV equation is not trivial. As for the problem of 
positivity, see Ref. [5]. 
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TABLE I: Comparison of the diffusion equation and the MCV equation 
III. MODIFIED MCV EQUATION 

So far, we have considered currents induced by gradients of n. However, when there exists 
a macroscopic velocity v, the total current is given by two contributions, 

Jtot=^V + J, (11) 

and the conserved density n should satisfy the equation of continuity with J tot instead of J 
itself. In this more general case, the MCV equation is modified |6]. 

Remember that the conserved density n, the velocity v and the thermodynamics force 
Vn are defined by the averaged quantities of particles which are contained inside a small but 
finite volume V*, commonly refered to as fluid cell [7]. Thus, Fick's law should be expressed 
for quantities per unit cell, 

[J(t)V*(t)) = -D[(Vn(t))V*(t)). (12) 

Note that V*(t) is a function of time because of the deformation induced by the macroscopic 
flow. In the case of Fick's law, however, V*(t) is a common factor and we can finally 
reproduce the ordinary result Q. On the other hand, Eq. ^ is modified as 

3(t)V*(t) = -J dsG(t- s)F(s)V*(s). (13) 

Because of the different time dependence of V*, the deformation of the fluid cell affects the 
evolution of the current J. The dynamics of the fluid cell is determined from a geometrical 
argument. Now we introduce three vectors to define the volume of the fluid cell, V* = 
(£i x £2) x £3- Because of the macroscopic velocity v, any vector is shifted to r — > r' = r + vdt. 
Finally we can derive the following equation for the deformation of the fluid cell, 

* V *(t) = (d t + v(t) ■ V) V*(t) = V*(t)V ■ v(t). (14) 
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Combining Eqs. (13) and (14), we obtain the modified MCV equation [6], 

TRj t J(t) + r R 3(t)(V ■ v(t)) + J(t) = —DF(t). (15) 

Interestingly enough, the final result does not depend on the volume of the fluid cell V*. 

The second term on the l.h.s. is a new term which does not exist in the MCV equation, 
and disappears when there is no macroscopic velocity. In the vanishing tr limit, this equation 
still reproduces Fick's law Q. 

The importance of the dynamics of volume elements is discussed also by Brenner [8] , where 
an additional velocity variable is introduced. In our case, however, we do not introduce any 
additional velocity. 

So far, we have discussed the modification of the diffusion equation by introducing a 
memory effect. There are, however, several different approaches to obtain modified diffusion 
equations. For example, there are attempts to obtain the MCV equation from the random 
walk model including memory effect, which is called the persistent random walk. So far, this 
approach reproduces the MCV equation only for the 1+1 dimensional case [9]. In our work, 
the memory function G has been assumed to have the exponential form. Other possibilities 
for the memory function are summarized in Ref. [TU]. Moreover, as is discussed in Ref. |llj . 
the problem of acausality may be solved by introducing a non-linear effect, van Kampen 
also discussed modifications of the diffusion equation, associated with the heat conduction 
of a photon gas [5]. The memory effect for phase separation is discussed in Ref. [12] 



IV. RELATIVISTIC FLUID DYNAMICS 



The problem of the violation of causality becomes more important in deriving relativistic 
fluid dynamics. In the following, we use the natural unit, where c = h = 1, and the metric 
^ = ^{1,-1,-1,-1}. 

In the derivation, we first choose gross variables which are necessary to extract the macro- 
scopic motion of many-body systems. If the chosen variables are not enough, the derived 
fluid dynamics will show unphysical behaviors, such as instability [EH E] and the divergent 
transport coefficients [3] . Unfortunately, there is no systematic procedure to collect the com- 
plete set of gross variables [15], but normally, conserved densities are chosen. In discussing 
phase transitions, order parameters are also one of the candidates of gross variables. Here 
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we do not consider phase transitions and conserved charges. Then the gross variables are the 
energy density e and fluid velocity ?P contained in the conserved energy-momentum tensor. 

In the idealized case, the the energy-momentum tensor is a function only of e and tP 
Then, by applying a Lorentz transformation and using the definition of the energy density 
and pressure P, we obtain T^ v — (e + P)u^u u — g^P. Note that P is calculated by the 
equation of state. Since T M " is conserved, we have 

<9 M T M " = 0. (16) 

This is the relativistic Euler equation. 

However, in general, T 7 *" cannot be expressed only by e and ip We represent this 
additional component by another second rank tensor IP". The most general T M " is, then, 
given by T^ v — (e + P)u^u u — g^ v P + IP* 7 . Conventionally, IP" is expressed using the trace 
part IT and traceless part 7r M " as IP" = 7r M " — (g^ u — u^u v )U [23]. Finally T M " is expressed 
as 

T M " = (e + P + U)u^u v - g^(P + IT) + tP", (17) 

and II and 7r^" are the bulk viscous pressure and the shear stress tensor, respectively, satis- 
fying the orthogonality condition u fJi TT flu = 0. 

The determination of these viscous terms is our next task. In the traditional Landau- 
Lifshitz theory [16] , these are induced instantaneously by the corresponding thermodynamic 
force, similarly to the diffusion equation, 

II = -£0 ?P" = 2rja^, (18) 

where ( and r\ are the bulk and shear viscosities, respectively. The thermodynamic forces 9 
and cP" are defined by 

6 = <VA (19) 
<P" = i (d^u" + 0V - ^{g" u - u»u u )8j = /\^ x5 d x u s . (20) 

When we use these definitions of the viscous terms, we obtain the relativistic NS equation. 
Because of the instantaneous production of the viscous terms, this equation contains the 
propagation with infinite speed, as is shown later. 

As was discussed for diffusion processes, the violation of causality is solved by introducing 
the memory effect previously discussed. Because of the existence of the fluid velocity P 1 , the 
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modified MCV equation should be applied. The relativistic representation of the modified 
MCV equation is 

T R u»d li J + T R Je = -DF. (21) 
Thus the viscous terms satisfying causality are given by 

nju^U + TaUe + U = -(9, (22) 
T n A^ x5 u a d a 7c xs + t v *TQ + if = 2tj<t^, (23) 

where Tn and are the relaxation times of II and tt^, respectively. Here the projection 
operator A^ A<5 is necessary to satisfy the orthogonality relation. In the following, we call 
this theory the causal dissipative relativistic (CDR) fluid dynamics. 



The second terms on the l.h.s. of Eqs. (22) and (|23[) comes from the deformation of the 



fluid cells. As is shown in Figs. 3-8 of Ref. j6], these terms are necessary to implement 
stable numerical calculations with ultra-relativistic initial conditions. 

Here, we simply assume that the thermodynamic forces of the modified MCV equation are 
the same as those of the relativistic NS equation. However, it may be possible to consider 
the higher order corrections to the thermodynamic forces, as is discussed in the Burnett 
equation [T7] . 



V. CAUSALITY AND STABILITY OF RELATIVISTIC FLUIDS 

When relativistic fluids are described by the relativistic NS equation, the fluid is Newto- 
nian because there is a proportional relation between tc' 11 ' and . On the other hand, the 
CDR fluid dynamics describes relativistic non-Newtonian fluids because tt^ is determined 
by solving the differential equation. 

Then, are relativistic fluids Newtonian or non-Newtonian ? To answer this question, we 
investigate the stability of these theories [TBI [Ti] . Let us introduce a perturbation ~ e iut-ikx 
around the hydrostatic equilibrium for the 1+1 dimensional system, 

e = e + 8e e iujt ~ lkx , U = 5U e iut - ikx , ~ = 5~ e iult - lkx , (24) 

where Eq = const and the fluid velocity is parameterized with S as / = (coshH, sinhS). 
Then the linearized CDR fluid dynamics is summarized as 

AX = 0, (25) 
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FIG. 3: The imaginary part of the dispersion relations for the relativistic NS theory (left panel) 
and the CDR fluid dynamics (right panel) at the rest frame |13j . The relativistic NS theory has 
two solutions, while the CDR fluid dynamics has three solutions. 



where 



X = (5e,5E,5U), 
I 



A 



\ 



iu -ik(e + P) 
— ikc 2 s iu(e + P) —ik 
ik( 1 + iuTu 



(26) 
(27) 



where P$ = P(so), and c s is the velocity of sound. In the following calculation, we consider 
a massless ideal gas, where Eq = 3Pq and c 2 s = dP/de = 1/3. Note that the result of the 
relativistic NS equation is reproduced by taking r n = 0. 

The dispersion relation is obtained from det A = 0. When rn ^ 0, we have one non- 
propagating mode and two propagating modes. From the propagating modes, the group 
velocity is calculated as 

Va = ^r~f- + M^y (28) 

One can see that the group velocity diverges in the vanishing rn limit. That is, the relativistic 
NS theory contains infinite velocity propagations. 

In order to study the stability, the imaginary part of the dispersion relations is shown in 
Fig. [3j We used C,/s = 0.1 and (C/ s )/( r n//3) = 1/6 where s is the entropy density and (3 is 
the inverse of temperature. In this parameter set, the group velocity is Vq = l/v^, which 
is a causal parameter set because the speed is slower the speed of light. The relativistic NS 
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FIG. 4: The imaginary part of the dispersion relations for the relativistic NS theory (left panel) 
and the CDR fluid dynamics (right panel) at the rest frame [13 . In this case, both theories have 
three solutions. One of the solution of the relativistic NS theory is negative. 

equation (left panel) has two solutions and the CDR fluid dynamics (right panel) has three 
solutions. All the imaginary parts are positive and both theories are stable. This result was 
already known in Ref. [T8] . 

If the theories are consistent with the relativistic kinematics, the nature of stability should 
not be changed by the Lorentz transform. To see this, we study the stability from a Lorentz 
boosted frame with a boost velocity of V = 0.1. As is shown in Fig. |4| one of the imaginary 
parts of the relativistic NS equation (left panel) becomes negative. On the other hand, the 
imaginary parts of the CDR fluid dynamics (right panel) are always positive. That is, the 
relativistic NS equation is inconsistent and inadequate as the theory to describe relativistic 
dynamics. 

These results suggest that the violation of causality and instability are correlated. To 
confirm this, we study the stability of the CDR fluid dynamics with an acausal parameter 
set, ((/s)/ (rn//3) = 1, where the group velocity exceeds the speed of light, vq = The 
result is given in Ref. [6]. Again, one of the imaginary part becomes negative in a Lorentz 
boosted frame, although all the imaginary parts are positive in the rest frame. That is, the 
violation of causality and instability is intimately related, and it is concluded that relativistic 
fluids are non-Newtonian. The results are summarized in table [IT} So far, we have discussed 
the bulk viscous pressure. The same result is obtained even for the shear viscous tensor [14]. 
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RNS CDR (acausal parameter) CDR(causal parameter) 
Rest frame stable stable stable 

Boosted frame unstable unstable stable 

TABLE II: Relation between causality and stability [13j 

VI. TRANSPORT COEFFICIENTS OF CDR FLUID DYNAMICS 

In fluid dynamics, transport coefficients are inputs which should be calculated from the 
underlying microscopic dynamics. In classical and non-relativistic NS fluids, it is known 
that the shear viscosity, for example, shows the following density dependence [26j[19j, 

V = Vo + ViP + W 2m P H • (29) 

The first term 770 can be calculated from two different approaches: the Chapman-Enskog 
expansion of the Boltzmann equation and the Green-Kubo-Nakano (GKN) formula. It is 
known that the both results are consistent. On the other hand, the higher order coefficients 
r]i and 772 are not calculated from the Boltzmann equation and we should use the GKN 
formula |27|. 

However, we cannot use the GKN formula to estimate the transport coefficients of the 
CDR fluid dynamics, because the GKN formula is derived by assuming the fluid to be 
Newtonian. Thus we have to derive a new formula to calculate the transport coefficients of 
the CDR fluid dynamics. 

The formula is derived by using the projection operator method [2U [22] ■ The results are 
summarized as 

V _ VGKN Tn 

f3(e + P)~ pzf d 3 x(f °*(x), f °»(0)) ' /3 

C (gkn ni 

0{e + P) ~ p j d 3 x (f o*( x ) ? f °*(0)) ' (3 

where " denotes operator, LT = Yl^=i T u /3— c 2 s T m and 5 A = A— Tr[p eq A] with the equilibrium 
density matrix p eq . Here t]gkn and (gkn are the shear and bulk viscosities of Newtonian 
fluids which are calculated using the GKN formula (more exactly, the Zubarev method). 
One can see that the new transport coefficients are still calculated from the GKN formula 
with the normalization factors given by the static correlation functions. 



, (30) 

/3 2 /rf 3 x(T^(x),T^(0)) 

(GKN /„jx 

/3 2 /rf 3 x(5n(x),5n(o))' 
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FIG. 5: The temperature dependences of r]/(r n (e + P)) (left panel) and £/(7n(e + P)) (right panel) 
for pions, m = 140 MeV [22]. The solid, dashed, and dotted lines represent the results of the 
microscopic formula, the new kinetic calculation, and the IS calculation, respectively. 

The same coefficients can be calculated from the Boltzmann equation with Grad's moment 
method [23, 24J. Now we compare the new formula with the Boltzmann equation. For this 
purpose, we calculate the quantities 7}/{r^{e + P)) and £/(Tn(e + P)), because these are 
independent of the choice of the collision term of the Boltzmann equation. The behaviors 
of r]/(r n (e + P)) and C/( r n(£ + P)) are shown on the left and right hand sides of Fig. |5j 
respectively. The results from the new formula are shown by the solid lines. Note that, as 



was discussed in Eq. (29), the Boltzmann equation can be consistent with the microscopic 



formula only in the dilute gas limit. Thus we calculate the ratios in the leading order 
perturbative approximation. 

So far, two different results for the same ratios are known from the Boltzmann equation; 
one is the well-known result obtained by Israel and Stewart (IS) [23], and the other the result 
obtained recently by Denicol et al. (DKR) [21] [28]. The former and latter are plotted by 
the dotted and dashed lines, respectively. The DKR results predict larger ratios than those 
of IS, but are still smaller than the results of the microscopic formula. 

This difference comes from the effect of quantum fluctuations. In order to incorporate the 
quantum effect in the Boltzmann equation, the collision term is modified. The ratios are, 
however, independent of the collision term and hence quantum corrections are not included. 
Thus, to compare the two results, the effect of quantum fluctuations should be neglected. 
Then we find that the solid lines agree with the dashed lines [22J. As is shown in Fig. [HJ the 
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effect of these fluctuations is quantitatively large and cannot be ignored even in the high 
temperature limit. 

VII. SUMMARY 

We discussed the infinite propagation speed of the diffusion equation and introduced 
the Maxwell-Cattaneo-Vernotte (MCV) equation to solve this problem. The drawback and 
advantage of the diffusion and MCV equations are summarized in table [Tj The MCV equation 
should be modified when there exists a macroscopic flow. 

By using this modified MCV equation, we derived a relativistic fluid-dynamical model 
called the causal dissipative relativistic (CDR) fluid dynamics. On the other hand, another 
relativistic fluid model was obtained by the relativistic generalization of the Navier-Stokes 
(NS) theory. To see which theory is more adequate, the stability of the theories was studied 
with linear analysis. Then we found that the violation of causality and instability are 
intimately related and the theory becomes unstable if it contains acausal propagations. The 
relativistic NS theory contains such propagation, and hence is unstable, while the CDR fluid 
dynamics is causal and stable. In this sense, all relativistic fluids must be non-Newtonian. 

Finally, we discussed the calculation of transport coefficients of the CDR fluid dynamics. 
Because of their non-Newtonian nature, the Green-Kubo-Nakano (GKN) formula is not 



applicable. The formulae for the CDR fluid dynamics are shown in Eqs. (30 ) and (31 ). These 
formulae are consistent with the results of the Boltzmann equation with Grad's moment 
method. 
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